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Using a data sample of 1.06 x 10 8 tp' events collected with the BESIII detector in 2009, the 
branching fractions of XcJ pnn~ and XcJ pnn~n° (J=0,l,2) are measurer!*]. The results for 
XcO pnir~ and Xc2 — > pnir~ are consistent with, but much more precise than those of previous 
measurements. The decays of \ci — > pnit~ and XcJ — > pnir~Tv° are observed for the first time. 

PACS numbers: 14.20. Gk, 13.75. Gx, 13.25. Gv 



1 Throughout the text, inclusion of charge conjugate modes is implied if not stated otherwise. 
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I. INTRODUCTION 



Exclusive heavy quarkonium decays constitute an important laboratory for investigating perturbative Quantum 



Chromodynamics (pQCD). Compared to J/if) and ip' , relatively little is known concerning \cj decays [1|. More 
experimental data on exclusive decays of P-wave charmonia is important for a better understanding of the nature 
°f XcJ states, as well as testing QCD based calculations. Although these states are not directly produced in e + e~ 

n 

collisions, they are produced copiously in radiative decays of 4>' , with branching fractions around 9% 1]. The large ip 
data sample taken with the Beijing Spectrometer (BESIII) located at the Beijing Electron-Positron Collider (BEPCH) 
provides an opportunity for a detailed study of XcJ decays. 

Previous studies indicate that the Color Octet Mechanism (COM) 2j], ccg — > 2(qq) could have large effects on the 
observed decay patterns of these P-wave charmonia states 3j-|7j]. To arrive at a comprehensive understanding about 
P-wave dynamics, both theoretical predictions employing the COM and new precise experimental measurements for 
XcJ many-body final states decays are required. 

Also, the decays of XcJ with baryons in the final states provide an excellent place to investigate the production 
and decay rates of excited nucleon TV* states, which are a very important source of information for understanding the 
internal structure of the nucleon 8j . 

n 

The decays XcJ — > pnn~ and XcJ — > pnn~n° were studied by the BESH experiment with 14 x 10 6 ip' events [9J. Due 
to limited statistics and the detector performance, only Xco ~~ pnir~ and Xc2 —> prm~ were observed and measured 
with large uncertainties. A series of three-body and four-body decays of XcJ, includingchannels with similar final 



states to ours such as ppn° and ppn°Tr , were measured by the CLEO-c Collaboration 
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ill- 



In this paper, we present a measurement of XcJ decaying into pnir and pnir ir°. The samples used in this analysis 



consist of 156.4 pb 1 of ip' data corresponding to (1.06±0.04) x 10 8 
42.6 pb -1 of continuum data taken at -y/s = 3.65 GeV/c 2 . 



12| events taken at yfs = 3.686 GeV/c 2 , and 



II. DETECTOR AND MONTE-CARLO SIMULATION 



BESHI 13( is a major upgrade of the BESH experiment at the BEPCH accelerator 



141 for studies of hadron 



spectroscopy as well as r-charm physics [15|. The design peak luminosity of the double-ring e + e collider, BEPCH, 



is 10 cm s at center-of-mass energy of 3.78 GeV. The BESIII detector with a geometrical acceptance of 93% 



of 47r, consists of the following main components: 1) a small cell, helium-based main draft chamber (MDC) with 43 
layers. The average single wire resolution is 135 /iin, and the momentum resolution for 1 GeV/c charged particles in 
a 1 T magnetic field is 0.5%; 2) an electromagnetic calorimeter (EMC) comprised of 6240 Csl (Tl) crystals arranged 
in a cylindrical shape (barrel) plus two end-caps. For 1.0 GeV/c photons, the energy resolution is 2.5% (5%) and 
the position resolution is 6 mm (9 mm) in the barrel (end-caps); 3) a Time-Of-Flight system (TOF) for particle 
identification (PID) composed of a barrel part constructed of two layers with 88 pieces of 5 cm thick, 2.4 m long 
plastic scintillators in each layer, and two end-caps with 48 fan-shaped, 5 cm thick, plastic scintillators in each end- 
cap. The time resolution is 80 ps (110 ps) in the barrel (end-caps), corresponding to a K/tt separation by more than 
la for momenta below about 1 GeV/c; 4) a muon chamber system (MUC) consisting of 1000 m 2 of Resistive Plate 
Chambers (RPC) arranged in 9 layers in the barrel and 8 layers in the end-caps and incorporated in the return iron 
yoke of the superconducting magnet. The position resolution is about 2 cm. 

The optimization of the event selection and the estimation of backgrounds are performed through Monte Carlo 
(MC) simulation. The GEANT4-based simulation software BOOST 16| includes the geometric and material description 



of the BESIII detectors and the detector response and digitization models, as well as the tracking of the detector 
running conditions and performance. The production of the ip 1 resonance is simulated by the MC event generator 



KKMC 



17 1, while the decays are generated by evtgen [18[ for known decay modes with branching fractions being set 



to world average values [1 



, and by lundcharm 



1 91 ] for the remaining unknown decays. 



III. EVENT SELECTION 



Charged-particle tracks in the polar angle range | cos6*| <0.93 are reconstructed from hits in the MDC. Only the 
tracks with the point of closest approach to the beamline within ±5 cm of the interaction point in the beam direction, 
and within 0.5 cm in the plane perpendicular to the beam are selected. The TOF and dE/dx information are used to 
form particle identification (PID) confidence levels for tt, K and p hypotheses. Each track is assigned to the particle 
type that corresponds to the hypothesis with the highest confidence level. Exactly one proton and one 7r _ , or one 
antiproton and one n + in the event are required in the analysis. 

Photon candidates are reconstructed by clustering the EMC crystal energies. The minimum energy is 25 MeV 
for barrel showers (|cos0| < 0.80) and 50 MeV for end-cap showers (0.86 < |cos6*| < 0.92). EMC cluster timing 
requirements are made to suppress electronic noise and energy deposits unrelated to the events. 
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A. Selection of ip' — >• "iXcJiXcJ — * pnir 

For the channel ip' — > JXcJ,Xcj — > pnir~ , at least one photon with energy greater than 80 MeV is required. To 
remove photons coming from interactions of charged particles or neutrons (antineutrons) in the detector, the angles 
between the photon and the antiproton, antineutron and other particles (pion, proton and neutron) are required to be 
greater than 30°, 20°, and 10°, respectively. A one-constraint (1C) kinematic fit is performed under the ip' — > jpnir~ 
hypothesis, where the neutron (antineutron) is treated as a missing particle. For events with more than one isolated 
photon candidate, the combination with the smallest x\c ^ s retained, and x\c 10 ^ s required. The distributions of 
the recoiling mass against 7p7r~ and jpTT + are shown in Fig. [TJ (a) and Fig. [TJ (b), respectively. Clear neutron and 
antineutron peaks are observed around 0.938 GeV/c 2 . 
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FIG. 1: The distributions of mass recoiling (a) against ypir~ in ip' — > ^pnn" and (b) against 7p7r + in ip' —¥ 7pn7r + . 
Dots with error bars are data, and the solid histograms are the sum of signal and backgrounds, where the 
backgrounds are estimated from the inclusive MC and the continuum data at -^/s = 3.65 GeV/c 2 . The dominant 
background contributions from ip' —¥ pmr - (pmr + ) and ip' — > pnir~ ir° (pnir + ir°) are shown as dashed and dotted 
lines, respectively. 



An antineutron can form a cluster in the EMC with very high probability due to annihilation in the detector. To 
further purify the events with antineutrons and more than one photon, a < 15° is required, where a is the angle 
between the expected antineutron direction and the nearest photon. Further, M pn -(p^+\ > 1.15 GeV/c 2 is required 
to remove background events with A or A. Finally, the transverse momentum for the proton or antiproton is required 
to be greater than 0.3 GeV/c due to the difference of the tracking efficiency at low transverse momentum between 
data and MC. 
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B. Selection of ip' — »■ 7Xcj, XcJ pnn 7T° 

For ip' — » 7X c j, XcJ — > pnir~Tr° , at least three isolated photons are required, where the photon isolation criteria are 
the same as those in ip' — > jXcj, XcJ ~^ pnir~. 7r° candidates are selected from any pair of photon candidates that can 
be kinematically fitted to the ir° mass and satisfy x 2 < 20. There must be at least one 7r° candidate. A 1C kinematic 
fit is performed to the ip' — > jjjpnn^ hypothesis constraining the mass of the missing particle to that of the neutron, 
where the three photons are a 7r° candidate together with another photon. Kinematic fits are carried out over all 
possible three-photon combinations. The kinematic fit with the smallest y 2 _ _ is retained, and y 2 _ _ < 10 is 
required. If there is more than one 7r° candidate among the three photons, the pair with the invariant mass closest 
to the 7r° mass is assigned to be from the ir° decay. To further suppress backgrounds with final states 77pmr~ and 
7777pn7r - , 1C kinematic fits are performed under 77pn7r _ and 7777pn7r~ hypotheses, 1 and x 7 fi3r - < X^-ypn^- an d 
X 2 - - < Y 2 - - are required. With the above selection criteria, the invariant mass of tt° and the distribution 

777P n7r /Vj 7777pn7r ^ ' 

of the recoiling mass against 7p7r~7r° in ip' — > 7pfi7r~7r are shown in Fig. [2] (a) and Fig. [5] (b), respectively. Similar 

distributions are obtained in the charge conjugate channel. 
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FIG. 2: The invariant mass distribution of (a) 77 from n° and (b) the distribution of recoiling mass against 
7P7T~7r° in ip' — > jpfiTv~Tv . Dots with error bars are data, and the solid histograms show the sum of signal and 
backgrounds, where the backgrounds are estimated from the inclusive MC and the continuum data at y/s = 3.65 
GeV/c 2 . The dominant background contributions from ip(2S) — > pn-K~ir° and ip(2S) —¥ n°n J/tl>,J/ip — > pnw~ 
are shown as dashed and dotted lines, respectively. 

To suppress the background from ip' — > 7r°7r° J/tp, for events with at least four photons, 7r°7r° combinations are 



1 The two-photon combinations are pairs of photons from 7T° candidates, and the four-photon combination are pairs of photons from n° 
candidates together with another two photons. If there are multiple combinations, the one with the smallest x 2 i s selected. 



formed with any four photons, the one with the smallest A = \J (m 7 i 7 2 — m n o) 2 + (m 7 3 7 4 — m n o) 2 is selected, and 



\M r 



M 



> 50 MeV/c is required, where M reco u is the mass recoiling against 7r°7r°. To purify the channel 



with antineutrons, the same requirement a < 15° as that for the channel tp' — ¥ ^XcJi XcJ ~^ pnn~ is applied, and 
the transverse momentum for proton or antiproton is required to be greater than 0.3 GeV/c in order to reduce the 
systematic uncertainty. 

With the above criteria, clear XcJ signals are observed in the invariant mass distributions of pnir~ (pnir^) in 
XcJ pnir~ (pnir + ) and of pfm~ ir° (pnir + ir°) in XcJ — > pfnr~ir Q (pnir + ir°), as shown in Figs.[3](a) and (b) and Figs.0] 
(a) and (b), respectively. 
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FIG. 3: Invariant mass distribution of (a) pnix~ for if)' — > r ypfnr~ events, and of (b) pnir + for the charge conjugate 
channel. Dots with error bars are data, and the filled histogram is the normalized non-resonant contribution 
estimated from continuum data at y/s = 3.65 GeV/c 2 . Resonant background, shown as the dashed histogram, is 
dominant and is estimated from MC. The sum of both background contributions is shown as the solid histogram. 



IV. BACKGROUND STUDIES 

Background is investigated using both exclusive and inclusive MC samples, as well as continuum data. For tp' — > 
IXc.Ji XcJ pnir~ , a number of individual channels from the inclusive MC sample have been found to contribute to the 
background. The following reactions have been found to be important: tp' —> pnir~ir° (sa 20% of the total background), 
tp' — > pwr~ (8%), tp' -> Tr°Tr°J/tp with J/tp -> pnir~ (1%), and tp' -> j\cj with XcJ -> pnir~Tr° (4%). For the first 



,t 



are used for normalization, and 



three background channels, the world average branching fractions listed in PDG 
for XcJ pnTT~TT°, we use our own measurement as described in Sect. IVII The background originating from non- 
resonant processes is estimated using a continuum data sample collected at a center-of-mass energy of 3.65 GeV/c 2 
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FIG. 4: Invariant mass distribution of (a) pfnr~iT for the if)' — > ^ypfnv~-K° channel, and of (b) pn7r + 7r° for the 
charge conjugate channel. Dots with error bars are data, and the filled histogram is the normalized non-resonant 
background contribution estimated from continuum data at *Js = 3.65 GeV/c 2 . Resonant background, shown 
as the dashed histogram, is dominant and is estimated from MC. The sum of both background contributions is 
shown as the solid histogram. 



after normalizing it to the integrated luminosity and the production cross section. The invariant mass distributions of 
pnn~ and of the charged conjugate state pnir + are shown in Figs. [3] (a) and (b), respectively. No peaking background 
is observed in the signal region. 

Also ip' — > JXcJtXcJ ~~ ^ pnTr~ir° contains a number of individual background channels according to studies with 
the inclusive MC sample. The dominant sources are the reactions ip' — > pnir~ir° (« 25% of the total background), 
ip' — » TT°7T°J/tp with J/tp — > phir~ (10%), and i/)' — > IXcj with Xc.j pnir~ (1%). Again, for the latter channel 
the branching fraction obtained from our own measurement is used for normalization (see Sect. I VII) ; for the first two 

n 

background channels the PDG [l| branching fractions are used. Additional background is studied using the inclusive 
MC sample, and the continuum data sample. The invariant mass distributions of phir~ir° and pnir + ir° are shown in 
Fig. @] (a) and Fig. 2] (b), respectively. Also in this case, no peaking background is observed in the signal region. It 
has been explicitly verified that the process "0' — > JXcj,Xcj "fJ/ip with J/tp — > pnir~ or J/tp — > pnir~ir° does not 
represent a significant source of background. 



V. INTERMEDIATE STATES 



We have searched for potential intermediate N* baryon resonances decaying into either pir~ or nir~ . Such a study 
is needed for correct modeling of Xc.i pnir~ and XcJ ~ * pnir~ir a in order to determine the efficiencies used in 
the calculation of the branching fractions of these decays, which is the main purpose of this work. Selecting Xco 
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events in the pnir~ decay mode, \M pfi7T ~ — M Xc0 \ < 45 MeV/c 2 , the efficiency corrected pir~ and nir~~ invariant mass 
distributions are shown in Fig. [5] (a) and Fig. [5] (b), respectively. Background contributions are also shown and are 
obtained from the sideband invariant mass region 45 MeV/c 2 < \M pfi7T ~ — M Xc0 \ < 75 MeV/c 2 . Fig. [5] (d) shows 
in addition the corresponding Dalitz plot M 2 (n7r~) vs. M 2 (p-K~). The structures of the N* states at around 1.4 
GeV/c 2 and 1.7 GeV/c 2 can be seen in both the pn~ and fvK~~ invariant mass spectra as well as the bands in the 
Dalitz plot. The pn invariant mass is shown in Fig.[S](c). A large enhancement in the pn threshold region is observed, 
which is also visible as a diagonal band along the upper right-band edge in the Dalitz plot in Fig. [5] (d). A similar 
thresho ld enha ncement has qualitatively been observed elsewhere, such as a pp threshold enhancement in B meson 
decays 



at BaBar 



27|. 



20j25j, V' decays 



26j , and in the shape of the timelike electromagnetic form factor of the proton measured 



The peak at around 2.2 GeV/c 2 in both the pir~ and mr - invariant mass spectra is partly due to the reflection of the 
threshold enhancement of pn. It might also be partly due to high mass N* states, such as TV* (2190) or TV* (2220). The 
same structures are observed in the charge conjugate mode of Xco —> pn-K + and the decays of x&,2 ~^ pnir~ . A partial 
wave analysis is necessary to obtain more information about the N* components and the threshold enhancement in 
the invariant mass distribution of pn. 
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FIG. 5: The invariant mass distributions with the efficiency correction of (a) pn~ , (b) mr~ , (c) pn and 
(d) Dalitz plot for Xco — > pnn~ events. The dots with error bars are data, and the histograms are for 
backgrounds obtained from sideband events. 



In X 



cJ 



pnir 7T°, various two-body and three-body invariant mass distributions for events within the \cj signal 
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region were investigated as well. No obvious N* state is observed. The distributions are very similar to those of phase 
space, except the 7r ± 7r° invariant mass spectra that show a significant signal (see Fig. [6]). 



250 



> 200 

CD 

o 150 

CM 



CD 

W 50 



8 







! i 

niixu 


1 + Data 
'{ — PHSPMC 


" /■ i 


T ' + TT1+L 

i > t i t t t i i i r*p*sui — 



2 0.4 0.6 0.8 1 1.2 1.4 1.6 

2\ 



M(7t-7t°) (GeV/c 2 ) 

FIG. 6: The invariant 7r _ 7r° mass distribution in Xc.j — > pn-K~n where a significant p signal is observed in the 
data. 



VI. SIGNAL EXTRACTION 

Signal yields are extracted using unbinned maximum likelihood fits to the observed pnir~ and the pnir~ir° invariant 
mass distributions. The following formula has been used for the fit: 

2 

Y^,BW{m\M i ;T i )®G{m;(T i ) + BG, (1) 

where BW(m;Mi; Ti) is the Breit-Wigner function for the natural lineshape of the XcJ resonance, BG represents 
the background shape and is described by a third order Chebychev polynomial, and G(m; cr,) is a modified Gaussian 
function parameterizing the instrumental mass resolution, which was used by ZEUS Collaboration in ref 
expressed by: 



28] and 



i * 1 >■ i+| ™_ ' 

G{jn>,<Ti) = -=^e-^ CT< . (2) 

V 27T<7j 

In the fit, the natural widths of the XcJ states are fixed to the PDG 1[ values, while their masses and corresponding 
instrumental resolutions are floated. For tjj' — > JXcJiXcJ ~^ pnn~ , this fit is performed in the mass region of 3.30 
GeV/c 2 < M{pn-K~) < 3.60 GeV/c 2 , for the process -0' ~> IXcJiXcj — > pnir~ir° the mass region of 3.30 GeV/c 2 < 
3.64 GeV/c 2 has been chosen. 

The fits to the pn-K~ and pmr + invariant mass distributions are shown in Fig. [7] (a) and [7] (b), respectively. 
The corresponding fits to the pfiir~ir and pnir + ir° invariant mass distributions are shown in Figs. |5] (a) and (b), 
respectively. 
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MC samples for signal events have been generated to obtain the relevant detection efficiencies. In the underlying 
event generators, for the decay ip' — > "fXcj an angular distribution proportional to 1+A cos 2 (f9) has been assumed, 
where 9 is the angle between the direction of the radiative photon and the positron beam, and A = 1, -1/3, 1/13 
for J = 0, 1, 2, respectively, in accordance with expectations of electric dipole (El) transitions. Since in the case 
of ip' — > "/XcJtXcJ ~^ pnn~ , structures have been observed in the j>7r - , n7r _ , and pn invariant mass spectra (see 
Sect. IV]) . the decays of Xcj into pnir~ are generated taking these structures and the polar angle distribution of the 
proton/neutron into account. The MC samples of Xcj — >• pnir~ir° used to determine the detection efficiencies were 
generated with a flat angular distribution, although a p^ intermediate state is seen. To estimate the systematic 
uncertainty associated with the missing p^ intermediate state, a MC production of Xcj pnp~ with the correct 
angular distribution of p ± — > ^it has been generated (see Sect. IVIlT|) . 
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The XcOi Xci and Xci MC samples are finally weighted with the amplitudes observed in data, and the same fitting 
process as that for data is performed to the mixed MC sample. The detection efficiencies of Xcj are calculated by 
e c j = N^j* /Nj!j n , where A^fJ* is the number of Xcj events extracted from the fit, and N^j n is the number of generated 
Xcj events. 



VII. BRANCHING FRACTIONS 



A. Branching fractions of Xci — > pnn 



The branching fractions of Xcj pnn are calculated according to: 

B{XcJ -> pmr - ) = — — Nsw r , (3) 

N. r x B{4>' -> 7x cJ ) x e cJ 

where N S i g is the number of signal events extracted from the fit to the invariant mass distribution, N^i is the number 
of tp' events, B(ip' —> JXcj) is the branching fraction of ip' —> jXcJ a s quoted in the PDG [l|, and e c j is the detection 
efficiency. The results are summarized in the left column in Table |TJ The same calculation for the charge conjugate 
channel is performed, and the results are summarized in the right column. 

TABLE I: The number of signal events N a i g , the detection efficiency e c j, and the branching fractions of Xcj — > pfin~ , where 
the errors are statistical only. 







Xc j — > pn-K 






XcJ — > pmr + 






XcO 


Xci 


Xc2 


XcO 


Xci 


Xc2 


N si g 


5150±102 


1412±58 


3309±79 


5808±121 


1625±73 


3732±89 


EcJ (%) 


38.6±0.2 


35.9±0.3 


39.2±0.2 


40.9±0.2 


40.7±0.3 


41.2±0.2 


b (icr 3 ) 


1.30±0.03 


0.40±0.02 


0.91±0.02 


1.38±0.03 


0.41±0.02 


0.98±0.02 



B. Branching fractions of Xcj — > pnn ty° 



Considering the branching fraction of ir° — > 77, the branching fraction of Xcj ~~ ^ P n7T 71-0 is calculated according to: 
B( XcJ - P nn-n°) = N ^ ixB ^^ lx ^l B{n ^ llU£c/ (4) 
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The results are summarized in the left column in Table HT1 The corresponding results for the charge conjugate channel 
are also shown in the right column. 

TABLE II: The number of signal events N S i g , the detection efficiency e c j and the branching fractions of Xcj — > pnn~-K° , where 
the errors are statistical only. 





XcJ 


— y pnir 7T° 




XcJ 


— > pnir + TT° 






XcO 


Xcl 


Xcl 


XcO 


Xcl 


Xc2 


N sig 


2480±85 


f082±52 


2128±62 


2757±94 


126f±60 


2352±69 


EcJ (%) 


10.4±0.1 


10.4±0.2 


9.8±0.1 


12.2±0.1 


12.3±0.2 


11.2±0.1 


B (f(T 3 ) 


2.36±0.08 


f.08±0.05 


2.38±0.07 


2.23±0.08 


f.06±0.05 


2.30±0.07 



VIII. ESTIMATION OF SYSTEMATIC UNCERTAINTIES 

Several sources of systematic uncertainties are considered in the measurement of the branching fractions. These 
include differences between data and the MC simulation for the tracking algorithm, the particle identification (PID), 
photon detection, the kinematic fit, the requirement on the angle a, ir° reconstruction, the fitting procedure, and 
the number of ip' events. Also possible imperfections in the description of intermediate resonances in the MC are 
considered. 

a. Tracking and PID The uncertainties from tracking efficiency and PID are investigated using an almost 



background-free control sample of J/tp — > ppir + n~ from (225.3 ± 2.8) x 10 8 J/4 1 decays 29(. The tracking efficiency 
is calculated with e = Nf u u/N a u, where Nf u u is the number of events with all final tracks reconstructed successfully 
and N a u is the number of events with three of them reconstructed while one track is missing. The PID efficiency is 
the ratio of the number of selected events with and without PID. Both efficiencies are studied for pions and protons 
(antiprotons) as a function of transverse momentum and cos#. The data - MC simulation difference for the tracking 
efficiency is estimated to be 1% per track. Therefore, a 2% uncertainty is taken for two-track events. For the PID 
efficiency, a 2% difference between data and MC is found for antiprotons, and 1% for any other charged particle. 
Therefore, 2% (3%) is taken as the systematic uncertainty for the final states including pir~ (pTr + ). 

b. Photon detection The uncertainty due to photon detection and photon conversion is 1% per photon. This 
value is determined from a study using clean control samples, such as J/ip — > p°n° and e + e~ — > 77. Therefore, a 
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1% uncertainty is taken for the tp' — » 7Xcj — > jpnir~ channel, while for the ip 1 — » 7X C J — > 7pn7r~7r° channel, with 3 
photons in the final state, a 3% uncertainty is taken. 

c. Kinematic fit The systematic uncertainty stemming from the 1C kinematic fit is investigated using J/tp — ¥ pfnT~ 
and J/tp —> pnir~ir° events, where n is treated as a missing particle with a mass of 0.938 GeV/c 2 . To obtain the 
systematic uncertainty associated with the fit, pure control samples of J/tp — > pnir~ and J/tp — > pnir^TT are selected, 
and the 1C kinematic fit is applied to both charge states, both for data and MC. The efficiency of the 1C kinematic 
fit is estimated calculating the ratio of the number of events with and without the kinematic fit. From the two charge 
conjugate states, the larger difference between data and MC is taken as the systematic uncertainty. We assign an 
uncertainty of 2.9% for the kinematic fit in the case of pn7r~(7r°) and 2.7% for pnir + (Tr°). 

d. Angle a Another source of systematic uncertainty is the requirement on the angle a 7fi < 15°. This uncertainty 
is studied using the J/tp — > pnir~ sample and an uncertainty of 1.8% is assigned for this item. 

e. ir° reconstruction The uncertainty from the ir° reconstruction is determined with a high purity and high statistics 
sample of J/tp — > 7r + 7r~7r . The tt° selection efficiency is determined by counting the number of 7r° candidates in 
the 7r + 7r - recoiling mass distribution with and without the standard it selection (1C kinematic fit). The data - MC 
simulation difference has been measured to be 0.7%, and it has been verified that there is no dependence of this value 
from the 7r° momentum and ir° polar angle. 

/. Uncertainty from intermediate states As mentioned above, in XcJ ~^ pnn~ , obvious N* intermediate states 
and a pn threshold enhancement are observed in the invariant mass spectra of p-n~ (mr - ) and pn (Fig. [5]). To 
account for these structures in the determination of the detection efficiency for \ c j — > pfm~, MC samples were 
produced including these structures at the generator level. To determine the systematic error associated with this 
procedure, an alternative method is used, where the efficiency, including the effect of these structures for XcJ pnit~ , 
is determined by re-weighting MC events generated according to phase space by the ratio of data and MC events 
in the two-dimensional distribution of pn~ versus ni:~ invariant mass. The difference in the detection efficiencies 
determined by the two methods is assigned as the systematic uncertainty. In \ c j pfnr~ir° : a strong signal is 
observed in the 7r ± 7r° invariant mass spectrum (Fig.[S]), while MC samples used to estimate the detection efficiencies 
are generated assuming phase space only. To study the effect of the intermediate state on the efficiencies, a sample 
°f XcJ pnp~ events is generated and analyzed. The efficiency difference between the two MC samples with and 
without a p intermediate state, {e pri7r ^o — e pnp )/e pn7T T T o , is taken as the systematic uncertainty. 
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g. Fitting procedure As described above, the yields of the \cj signal events are derived from fits to the invariant 
mass spectra of pnir~ and pnir~ir°. To evaluate the systematic uncertainty associated with the fitting procedure, we 
have studied the following aspects: (i) Fitting range: In the nominal fit, the mass spectra of pfnr~ and pfnr^ir are 
fitted from 3.30 GeV/c 2 to 3.60 GeV/c 2 and from 3.30 GeV/c 2 to 3.64 GeV/c 2 , respectively. We have changed these 
intervals to 3.32—3.60 GeV/c 2 in the case of pfnr~ and 3.32—3.62 GeV/c 2 for pwr _ 7r° events. The differences in the 
finally obtained branching fractions by changing the fit intervals, are taken as the systematic uncertainties associated 
with the fit intervals, (ii) Signal lineshape: The partial width for an El/Ml radiative transition is proportional to 
the cube of the radiative photon energy (E^), which leads to a diverging tail in lower mass region. Two damping 



31] collaborations and have been used in addition 



functions have been proposed by the KEDR |30( and the CLEO 
to the standard approach in the formula describing the fit to the signal lineshape. Differences with respect to the 
fit not taking into account this damping factor have been observed, and the greater of the two differences is taken 
as the systematic uncertainty associated with the signal lineshape. (iii) Mass resolution parameterization: A single 
Gaussian formula instead of the modified Gaussian formula (Formula 2) has been used to describe the instrumental 
mass resolution. The resulting differences for the final branching fractions are taken as the systematic uncertainties, 
(iv) Mass resolution: Studies have shown that the XcJ mass resolutions, as simulated by MC, are underestimated. 
To evaluate the systematic effects associated with this aspect, the invariant masses of pnir~ and phir~Tr° in the MC 
samples are smeared with a Gaussian function, where the width of this Gaussian depends on the invariant mass as 
well as on the channel. The same fitting processes as in the nominal cases are performed on the smeared mass spectra 
of pnir~ and pnir~ir°, and the detection efficiencies are recalculated. The efficiency difference between the smeared 
and unsmeared case is taken as the systematic uncertainty, (v) Background shape: To estimate the uncertainties due 
to the background parameterizations, a second order instead of a third order Chebychev polynomial is applied in the 
fitting. Again, the difference between the two cases is used as an estimate of the systematic uncertainty. 

h. Other systematic uncertainties The number of ip' events is determined from an inclusive analysis of ip 1 hadronic 



events and an uncertainty of 4% [12J is associated to it. The uncertainties due to the branching fractions of ip' — > ^fXcj 
and 7T° — > 77 are taken from the PDG [l| . A small uncertainty due to the statistical error of the efficiencies is also 
considered. 

In Table IIIII a summary of all contributions to the systematic error is shown. In each case, the total systematic 
uncertainty is obtained by adding the individual contributions in quadrature. 
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TABLE III: Summary of systematic errors (in %) for the branching fraction measurements of %cj — > pnir and \cj P n7T 71-0 • 





XcJ 


— > pnir 


XcJ 


— >• pnn + 


XcJ 


— > pTLTV 


n 
7T U 


XcJ ~ 


-> pnn 


4- 




XcO 


Xcl 


Xc2 


XcO 


Xcl 




XcO 


Xcl 


Xc2 


XcO 


Xcl 


Xc2 


MDC tracking 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


PID 


2.0 


2.0 


2.0 


3.0 


3.0 


3.0 


2.0 


2.0 


2.0 


3.0 


3.0 


3.0 


Photon detection 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


3.0 


3.0 


3.0 


3.0 


3.0 


3.0 


Kinematic fit 


2.9 


2.9 


2.9 


2.7 


2.7 


2.7 


2.9 


2.9 


2.9 


2.7 


2.7 


2.7 


Q 7 fi < 15° 


1.8 


1.8 


1.8 




- 


- 


1.8 


1.8 


1.8 








7r reconstruction 




- 


- 




- 


- 


0.7 


0.7 


0.7 


0.7 


0.7 


0.7 


Intermediate states 


5.3 


6.2 


8.2 


4.0 


1.2 


5.0 


3.0 


1.7 


0.7 


1.5 


1.8 


2.3 


Fit range 


0.3 


0.1 


0.1 


0.6 


0.2 


0.1 


0.0 


1.6 


1.1 


0.4 


2.4 


0.8 


Signal lineshape 


1.4 


2.4 


1.1 


0.9 


5.4 


1.4 


0.8 


2.5 


1.1 


1.5 


1.7 


0.8 


Resolution para. 


1.9 


4.6 


3.0 


3.1 


5.3 


3.1 


2.1 


5.5 


1.7 


0.4 


4.9 


2.2 


Resolution diff. 


U.o 


1.4 


1.3 


U.O 


0.2 


0.7 


1 o 

j-.y 


1 G 

i.y 


i n 


U.o 




1 Q 


Background shape 


1.7 


5.8 


1.2 


0.3 


1.1 


0.2 


1.3 


1.6 


1.1 


1.1 


2.4 


0.9 


N(i//) 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


B(ip' -> 7Xcj) 


3.2 


4.3 


3.9 


3.2 


4.3 


3.9 


3.2 


4.3 


3.9 


3.2 


4.3 


3.9 


B(ir° -> 77) 














0.03 


0.03 


0.03 


0.03 


0.03 


0.03 


MC statistics 


0.2 


0.3 


0.2 


0.2 


0.3 


0.2 


0.1 


0.2 


0.1 


0.1 


0.2 


0.1 


Total 


9.1 


12.5 


11.5 


8.6 


10.8 


9.5 


8.6 


10.6 


8.3 


7.9 


10.6 


8.8 



IX. SUMMARY 

Based on a data sample of 1.06x10 s ip' events collected with the BESIII detector, the branching fractions of 
XcJ pnir~ and \ c j — > pnTr~ir° are measured for J = 0,1,2. The results are summarized in Table |lVl where for 
each branching fraction the first error is statistical and the second systematic. The product branching fractions of 
B(ijj' — > 7Xcj) x B(\cj pnir~) and B(tp' — > jXcj) x B{XcJ ~ > pnir~ir°) are also summarized in Table Ivl For 
XcO pn-n~ and \ci — > pfnr~ , the results are consistent with the world average values within one standard deviation, 
while the precision is improved significantly. For the other XcJ decay modes, the branching fractions are measured 
for the first time. A comparison of individual branching fraction shows good agreement between charge conjugate 
channels. The measurements improve the existing knowledge of the XcJ states and may provide further insight into 
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their decay mechanisms. Based on the results of this work, detailed studies concerning the intermediate states may 
follow in the future. 

TABLE IV: Summary of branching fractions for XcJ ~* pnir~ and XcJ ~> pniT^ir . The first errors are statistical, and the 
second ones are systematic. 

XcO Xcl Xc2 

B{xc.J ~>pniv~) (1(T 3 ) 1.30±0.03±0.12 0.40±0.02±0.05 0.9I±0.02±0.I0 

B(Xc.j -> pnw+) (1CT 3 ) 1.38±0.03±0.12 0.41±0.02±0.04 0.98±0.02±0.09 

B{xc.J ^pn7r~7r°)(10 -3 ) 2.36±0.08±0.20 1.08±0.05±0.12 2.38±0.07±0.20 

BixcJ pn7r+7r°)(10" 3 ) 2.23±0.08±0.I8 1.06±0.05±0.11 2.30±0.07±0.20 

B{xc.j ~>pn-K~) (ICT 3 ) (PDG [1]) 1.14±0.31 - 1.10±0.40 



TABLE V: Summary of the product branching fractions of ip' — > JXcJ, XcJ P n7V and ip' ~~ > 7XcJ, XcJ ~~ > P n7T 71-0 • The 
first errors are statistical, and the second ones are systematic. 

XcO Xcl Xc2 

B{ip' iXcj) x B(XcJ -> pn-K~) (1CT 4 ) 1.26±0.02±0.11 0.37±0.02±0.04 0.80±0.02±0.09 

B(ip' -» iXcj) x B(xcj -»■ pn7r+) (1CT 4 ) 1.34±0.03±0.11 0.38±0.02±0.04 0.85±0.02±0.07 

6(V>' -»■ 7Xcj) x BixcJ pn-K-n Q )(lQ- 4 ) 2.29±0.08±0.I8 1 .00±0.05±0.I0 2.07±0.06±0.f5 

B{ip' -> iXcj) x B(xcj -> pn7r + 7r°)(f0" 4 ) 2.16±0.07±0.1 6 0.98±0.05±0.f0 2.0f±0.06±0.I6 
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